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Abstract 

We give a pedagogical introduction into an old, but unfortunately not very well-known for- 
mulation of GR in terms of self-dual two-forms due to Plebanski. Our presentation is rather ex- 
plicit in that we show how the familiar textbook solutions: Schwarzschild, Volkoff-Oppenheimer, 
as well as those describing the Newtonian limit, graviton and homogeneous isotropic Universe 
can be obtained within this formalism. Our description shows how Plebanski formulation gives 
quite an economical alternative to the usual metric and frame-based schemes for deriving Ein- 
stein equations. 
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1 Plebanski formulation of general relativity 



The aim of this short paper is to give a description of Plebanski formulation of general 
relativity (GR) in a version that we found most suited for practical computations. Our pre- 
sentation is very explicit, in that the standard textbook solutions of GR are obtained. As we 
shall see, given an antsatz for the metric, Plebanski formulation produces Einstein equations 
even more quickly than the already efficient tetrad method. In our opinion, the efficiency and 
beauty of this formulation may warrant its inclusion in general relativity textbooks. 

Our convention for the signature is (— , +, +, +). We start with a description of the Plebanski 
version of Einstein equations in the usual tensor notations. 



1.1 Einstein condition and the Hodge operator 

Given a spacetime metric g^i, the condition that this metric is Einstein reads: R^j^u ~ Qf^^, 
where R^^, := Rff is the Ricci tensor, and the fact that the proportionality coefficient in this 
condition must be a constant is implied by the (differential) Bianchi identity V^G^y = 0, where 
is the metric-compatible derivative operator and G^jy = R^^y — {l/2)gp_yR is the Einstein 
tensor. As usual, the quantity R is the Ricci scalar R := R'^, and all indices are raised and 
lowered with the metric. The Plebanski formulation of GR is based on the following simple 
and well-known reformulation of the Einstein condition in terms of the Hodge operator. Thus, 
let us introduce the operation of Hodge dual that acts on bivectors (anti-symmetric rank two 
tensors) Ap^i,: 

A * A — -t^ P"^ A (^^ 

where the quantity e^j^po- is the volume 4-form for the metric g^^p. The following elementary 
properties of the Hodge operator are easily verified: its square is minus one and it is invariant 
under conformal transformations of the metric g ^^g- 

Given the Riemann curvature tensor Rp,upa one can apply the Hodge operator to either the 
first or the second pair of indices: 



1 „/,/„ 1 

2' 



J^pupa ■— r)^pu J^p'u'pa-, ^pupa ■~ ^ -'l-Ati^p' ct' t pa- ) 



It is a straightforward computation to check that the Einstein condition (together with the first 
Bianchi identity Rp,[upa] = 0) is equivalent to the condition that the left and right Hodge duals 
of the Riemann tensor coincide: 

Rpy ~ 9pii ^ Rpvpa — Rpupa- (3) 
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Indeed, since {*Y = —1 the last condition is equivalent to R/j^upa = ~*R*^upa- Contracting a 
pair of indices on the left-hand-side to produce the Ricci curvature, and expanding the product 
of two epsilons on the right-hand-side in terms of products of the metric tensor, one finds that 
this condition gives R^^, = {l/4)g^yR, which is the Einstein condition. 

To arrive to Plebanski formulation we need a slight additional reformulation of the condition 
([3]). Thus, let us use the Hodge operator ([T]) to introduce projectors on the spaces of the so- 
called self- and anti-self-dual bivectors. These are bivectors that remain essentially unchanged 
under the Hodge duality: 

*A^"' = iA""" ^ (self - dual), "A^" = -lA^'' ^ (anti - self - dual). (4) 

The corresponding projectors are: 



where we have introduced the identity operator in the space of bivectors: 

/ /"^ := <5l/'<5'^] = i i^P^" - ^'^5^) . (6) 



It is then easy to show that the condition ([3]) can be rewritten using the above projectors 
as follows: 

(Ricci) ~ (metric) p- RP+ = 0, (7) 

where the anti-self-dual projector is applied on the left and the self-dual projector is applied 
on the right. Another convenient way to state this is to say that the Einstein condition is 
equivalent to the statement that the self-dual part of the Riemann curvature with respect to 
the second pair of indices is self-dual with respect to the first pair as well. In the mathematics 
literature this is known as the At iyah- Hit chin- Singer theorem [2]. A proof of equivalence of ([7j) 
to (IHl) is an elementary exercise using the definition of the projector operators. We can now 
use the observation ([7]) as the basis of a new formulation of general relativity. 



1.2 Plebanski formulation: pure gravity 

Let us present Plebanski formulation as a concrete recipe for obtaining Einstein equations 
starting from a metric. We will then explain why this recipe is equivalent to (JTj). 

The first step of this formulation is to find a tetrad for the spacetime metric. Thus, one 
represents the metric as: ds"^ = 9^ ® d^^ij-, where 9^,1 = 0,1,2,3 are the tetrad one-forms 
and T]jj is the Minkowski metric. The next step is to construct two-forms 6^ A 6"^ and then 
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take their self-dual parts with respect to the indices IJ . More concretely, one introduces an 
arbitrary time plus space split of the internal indices / = (0,?),? = 1,2,3, and considers the 
following triple of two-forms: 

2 ^ ' 

As is not hard to check, the two-forms satisfy: 

i 



AE^ = S'^^/^d^x, S^AS^ = 0, (9) 
where the anti-self-dual forms are given by: 

:= A 9' + A 6^. (10) 

After the two-forms are written down, the next step is to find a connection that is 
"compatible" with the triple S', i.e. the connection that satisfies: 

dS* + e^^'^A^' A S'^ = 0. (11) 

This is an algebraic equation for the components of the connection A^ that has a unique solution. 
It is not hard to write down an explicit expression for A* in terms of the derivatives c/S*, but in 
practice it is easier to solve the equation ffTTj) for each problem at hand by expanding it in the 
coordinate three- forms. We shall denote the E-compatible connection by Aj]. It is not hard to 
verify that Ay. is just the self-dual part of the tetrad-compatible dO^ + cjj A 6'"^ = connection 
uj^'^ , i.e., A^ = iio^"^ — {\j1)t^^^uj^^ . This fact will be important below for our explanation of the 
Plebanski method from the point of view of (JTj), but it is not essential if one just want to follow 
the method. 

The final step is to compute the curvature 

F' = dA' + ^e'^f'A^ A A'' (12) 

of the connection A^- Curvature is a two- form, and can be split in the basis of self-dual and 
anti-self-dual two-forms. Thus, one can always write: 

f\Ay) = f'^j:^ + F^s^'. (13) 

The matrices F^^ , F^^ are the main quantities of interest for the Einstein equations can be 
written quite easily as the following ten conditions: 

Tr(F) = -A, F'^ = 0, (14) 

where A is the cosmological constant. As a bonus of this method, one not only obtains Ein- 
stein equations, but also automatically gets an expression for the part of the curvature not 
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constrained by Einstein equations - the Weyl curvature. Thus, the self-dual part \1'*-' of the 
Weyl curvature tensor is given by the tracefree part of the matrix F^^: ^"^^ = {F'^^)tf. 

The described method of obtaining Einstein equations is quite efficient for practical com- 
putations. The steps one has to take to arrive at (fT^ are similar to the steps one makes in 
the tetrad formulation. The advantage of working with two-forms instead of tetrads 9^ is 
that one has half equations to deal with at intermediate steps, at the expense of all quantities 
becoming complex. We shall see the power of this method below when we use it to obtain the 
standard solutions of GR. 

An explanation of this method from first principles, i.e. from equation ([7]) is as follows. As 
we have already said, the S-compatible connection turns out to be just the self-dual part 
of the tetrad-compatible spin connection uj^'^ . This means that the curvature F^{Aj:,) is just 
the self-dual part of the curvature tensor R^-^{lj). Now the second equation in f|T^ just says 
that the curvature F^Aj:) is self-dual as a two-form. However, because of the noted relation 
between and R^'^ this is equivalent to the condition that the self-dual part of R^"^ is self-dual 
as a two- form, which is just the condition (I7j). The first equation in (fT4l) can be shown to be a 
consequence of the Bianchi identity DaF^ = 0, analogous to how the proportionality coefficient 
in the Einstein condition R^^, ~ g^,^ is shown to be a constant related to the cosmological 
constant. For more discussion on abstract aspects of Plebanski formulation of GR the reader 
may consult [3]. 

1.3 Coupling to matter 

In the previous subsection we have described the vacuum theory. We now need to extend it to 
the case when a non-zero right-hand-side of vacuum equations is present. 

The coupling to matter is described in Plebanski formulation in complete analogy to the 
usual metric-based scheme. The main step is to split the stress-energy tensor of matter T^^ 
into its trace T = and the tracefree T^i, = T^i, — {l/4:)g^^T parts, and then form out of the 
tracefree part a 3 x 3 (complex) matrix 

yii _ j^p^i^^giM,.^ (15) 

where, as before S*,S* are the self- and anti-self-dual forms ([H]) and ( ITOl) . The non- vacuum 
Einstein equations f|T^ then take the following simple form: 

Tr(F) = -A - 2nGT, F'^ = -2TiGT^ . (16) 

For the ideal fiuid, which is what matters for most practical applications, we have T^i, = 
(p + P)u^Ui, + Pg^ui where p, P are the energy and pressure densities correspondingly, and 
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is the 4- velocity vector. The trace of the stress-energy tensor is given by T = (3P — p), the 
tracefree part T^,^ = (p + P){u^Ui, + {l/'i)g^y) and so we get: 



where u'' is just the spatial component of the normalized u^u^ = — 1 velocity 4- vector: 

u, = J-^ iidt)^ + u\dx'),). (18) 

In these two formulas \u\'^ = u'^u\ with, as usual, a sum over the repeated index i implied. 
Other types of matter can be described similarly. 

We are now ready to study some simple solutions of GR using the formulation described. 
We start with the graviton solution of the linearized theory. We shall set the cosmological 
constant A to zero from now on. 



2 Graviton solution 



The graviton is a solution of linearized around Minkowski spacetime vacuum field equations. 
The background two-forms are given by: 

= idt A dx' - ^e'^'^dx^ A dx''. (19) 
In Plebanski formulation the graviton is described as a perturbation: 

= + 5S\ (20) 
The perturbation two-forms can then be decomposed into the background two-forms: 

6j:' = b''j:l + b'^f:i. (21) 

In terms of the matrices b^\b^^ the reality conditions iQ become 

b'^ = b6'^ + W = {W)\ (22) 

where 6 is a real scalar and uj"^^ is an arbitrary (complex) anti-symmetric matrix. The second 
equation implies that the matrix of quantities W is real. 

The representation (1211) of a general perturbation immediately provides a convenient clas- 
sification of perturbations. Indeed, the matrix can be split into its symmetric traceless, 
anti- symmetric, and the trace parts, while the matrix uj"^^ is already anti-symmetric. Then the 
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transverse part of an symmetric tracefree matrix describes a tensor mode, while the transverse 
part of an anti-symmetric matrix describes a vector one. For example, perturbations of the 
tensor type that correspond to the gravitons are given by: 

StensT.' = (23) 

where h'^^ is a symmetric, traceless, transverse h^l = three by three matrix. 

Perturbations of the scalar and vector types can also be considered but the modes they 
describe are non-propagating. Further, some of these modes are gauge. For example, the 
modes described by a;*-' part of the perturbation are just the gauge modes corresponding to 
availability of (complexified) SO (3) transformations in this framework. Other modes are gauge 
for the infinitesimal diffeomorphisms, which in this framework are described by: 

^^E^ = C^E' = di^E\ (24) 

where is the Lie derivative along an infinitesimal vector field ^ and is the operation of 
the interior product a vector field with a form. It is quite easy to compute and see which 
modes can be set to zero by the diffeomorphisms. One finds that, for example, the gauge- fixed 
perturbation of the scalar type is described by: 

6scai^' = (j)^' + (25) 

where 0, ip are related to the usual Newtonian potentials $, \I' in an elementary way. We shall 
further consider scalar perturbations below when we study the Newtonian limit. 

Having determined the form of the perturbation that we would like to study, it is elemen- 
tary to obtain the linearized Einstein equations. Thus, the first step is to find the linearized 
connection such that: + e^^^SA^ A S'^ = 0. This is a simple exercise in algebra with 
the result being: 

5A' = {-h'^e^^^ + i{h'^)') dx^. (26) 

Here the prime denotes the time derivative. The corresponding linearized curvature is also easy 
to compute. We decompose the answer into the self- and ant i- self- dual background forms: 

d5A' = ^{h'^)"{^^ + S^) + ^A0{J:^ - S^') + i(/iy)V'='S^'. (27) 

Now setting to zero the anti-self-dual terms on the right hand-side gives the Einstein equations: 

{h'^y - Ah'^ = Uh'^ = 0, (28) 

whose solutions are plane waves - gravitons. As a bonus, we also get an expression for the 
matrix of the components of the self-dual part of the Weyl curvature of the graviton: 

^ij = 1 ((/^iiy ^_ ^j^ij _^ 2i{h'^)'e^''^) , (29) 
which, as is easy to check, is symmetric and traceless. 
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3 Schwarzschild solution 



In this case one is still interested in vacuum field equations. We start from the usual expression 
for a spherically-symmetric metric: 

ds'' = -f{r)de + g\r)dr^ + r^d^\ (30) 

where dVL^ is the usual metric on the unit sphere. A tetrad is given by: 

e' = f{r)dt, e' = g{r)dr, = rdO, e'^ = rsm{e)d(j). (31) 

A convenient set of self-dual combinations ([8]) is then: 

= ie*Ae"-e^Ae^ = ie* A - e"^ A e^ = ze* A e"^ - e'' A (32) 

We now have to find the associated su(2) connection A\ This involves some guesswork, 
with the answer being: 

.1 'if , , , .9 sm( 9)dd) .■^ dO , . 

= -^dt + cos(e)d(j), = = —, 33 

9 9 9 

where, as usual, prime denotes the derivative with respect to the r coordinate. It is not hard 
to verify that the above connection solves (ITTll . 



In practice finding the connection is the most time consuming task. After this is done, it 
only remains to compute the curvature. This is a simple exercise in differentiation. One gets: 

= dA^ + A^ AA^ = - (^-^^ dtAdr-(^l- sm{e)de A 

F^ = dA^ + A^AA^ = -^ sm( e)d(f) A dr - %dt A dO, (34) 

9 9 

= dA^ + A^AA'^ = -^dr A dO - % sm{0)dt A dcf). 

9 9 

In the final step of extracting Einstein equations one simply has to express the coordinate 
two- forms appearing in (IMl) in terms of the self- and ant i- self- dual forms = — (S*)*. Thus, 
our final expression for the curvature components is: 



2/5- \ 9 J 2r2 \^ g 



= - S^) + t(^' + ) ' (35) 

^9 r \g f 



(S^-S^) + ^(S^ + S^) . 



Sfi-V \g f 



7 



After one has an expression for the curvature in terms of the basis two- forms S% one can 
immediately read off Einstein equations. First, one should equate to zero all the anti-self-dual 
components of the curvature. This immediately gives (from F'^,F^): 

f q' 

^ + - = 0, (36) 

whose solution is the familiar / = 1/ g. The equation coming from the anti-self- dual part of 
reads: 



2 I 1 - 71 )= 0- (37) 



fg^g) ^ r'^^ g 

The other equation, namely the condition that the trace part of the self-dual part of the 
curvature is zero reads, with (136|) . (1371) taken into account, after some simple rewriting: 

'^-^y (38) 



g^ J V g^ / 

which immediately gives the famous: 

4 = 1--' (39) 

which also solves fl37|) . as is not hard to see. The sign of the integration constant here is chosen 
so that the spacetime obtained is that of positive total gravitational mass. The integration 
constant itself is chosen so that r = r+ is the place where g~'^{r) vanishes. 

We would like to emphasize how much more thought-economizing the above derivation was 
as compared to the usual (e.g. tetrad-based) derivation given in all the textbooks. Even the 
steps leading to the computation of the curvature involve just the half of operations needed 
in the tetrad scheme (in the above spherically-symmetric case the economy is probably not 
as good, as a large number of coefficients is anyway zero). However, what really is most 
economical about Plebanski formulation is that one can immediately read off the equations 
from the expression for the curvature, unlike in the usual tetrad-based scheme. Indeed, in 
the tetrad scheme one first has to form the Ricci tensor, and only after that write Einstein 
equations. In this last step one has to raise and lower indices, which introduces some tricky 
minus signs, making it a non-trivial exercise to get the right final expressions. In the case of 
Plebanski formulation this last step is absent altogether, the manipulations leading to Einstein 
equations are as algorithmic as all the previous ones, making it much harder to make a mistake. 



4 Volkoff-Oppenheimer solution 

In this section we switch on the stress-energy of matter. We have computed the curvature 
components in the previous section. The equations that one obtains from the anti-self- dual 
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components of the curvature are: 



^' + 4 ) =27rG(p + P), (40) 



(41, 

The equation one obtains from the trace of the self-dual part is given by: 

Let us now take twice the equation (HOj) summed with fHTj) . and subtract the result from the 
equation (H2l) . We get, after some simple rewriting: 

^ ^1 - ^1 +\(l-^]= 8nGp, (43) 



which is the non-zero p generalization of the equation fl38l) . It immediately gives: 

2m r 
g~^{r) = 1 , m(r) = / AnGpr"^. (44) 



(45) 



r 

The function / can then be found from e.g. equation ( HOi) . Substituting (H4l) we get: 

/' _ m + AnGPr^ 
f r(r — 2m) ' 



integrating which we get /. This ODE should be solved with the "boundary" condition P{R) 
1 — 2M/i?, where M, R are the total mass and radius of the spherical object in question. 



5 Newtonian limit 

Let us now consider the case of small (static) perturbations around the Minkowski spacetime 
background (IT^ described by two scalar potentials 0, ip: 

5B' = + (46) 

Here we assume the functions (f), ip to be only those of the spatial coordinates (time inde- 
pendent). It can be easily shown that all other scalar perturbations can be set to zero by a 
diffeomorphism and an SO (3) rotation. It is easy to check that fH6l) corresponds to the following 
perturbation of the metric: 

rfs^ = -(1 + 2$)dt^ + (1 - 2^) dxl (47) 
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with 



^ (p + S-ip ^ - (j) 



We now have to solve the hnearized "compatibihty" equation and obtain the hnearized 
connection. We can solve for the (f) and ip parts of the perturbation separately. We get: 

4 = ^'PA - l^'"'<P,kdx\ = ^-^i;,dt + ^e*^V,fcC^a;^- (49) 

The corresponding linearized curvatures are: 

dA'^ = -^A0E^ + ^(-20,, + 5'^A<P)f:^, (50) 



Let us now analyze the Plebanski equations f|T6l) with zero pressure P = and velocity 
= 0. First, the anti-self-dual part on the right-hand-side is diagonal, which requires of the 
off-diagonal components to vanish: 

{{<P + ij),,)tf = 0, (51) 

where tf denotes the trace-free part. This implies + = 0, or \1' = $. Both of the diagonal 
components then give the Laplace equation, which in terms of the potentials (jlHl) takes the 
familiar form: 

A$ = 47iGp. (52) 
We also get for free an expression for the components of the Weyl tensor: 

^ij = _ -S'^A<^, (53) 
3 

which, we note, is real, as is typical of static backgrounds. 



6 Homogeneous isotropic Universe 

Using Plebanski formalism it is also quite easy to obtain the Friedman equations describing the 
evolution of a homogeneous isotropic Universe. The corresponding metric is given by: 



ds^ = a^ri) (^-dri^ + J^idx'Y^ 



(54) 
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where rj is the conformal time. The corresponding self- and anti-self-dual two-forms are given 
by: 

= a\idr] A dx' - ^e'^^dx' A dx^), = a^i.id'q A dx' + ]^e'^''dx^ A dx^). (55) 

The connection A* compatible with the set of self-dual two- forms S*, i.e. satisfying (|TT1) . is 
given by: 

A' = mdx\ (56) 

where we have introduced: 

a 

and the prime denotes the derivative with respect to the conformal time 
dA"- + {l/2)e''^^A^ A A^ of the connection ( l56l) is given by: 

F' = -^{w + 7^=')s' + -^{W - i-e)t.\ 

From fl58l) . f|T6|) we immediately get Einstein equations: 

^in' + n') = ^{p- 3P), ^{W - re) = -2nG{p + p). (59) 

These are more customarily known as their linear combinations: 

= 2n' + n' = -sirGaY (eo) 

3 

7 Discussion 

We would like to conclude this short paper by emphasizing once more how computationally- 
efficient the Plebanski formulation of general relativity is as compared to other standard text- 
book methods such as that based on tetrads. For this reason its inclusion in GR textbooks 
may be warranted, also in view of its conceptual simplicity and beauty. 

What we have not discussed in this paper is if Plebanski formulation can be not just a tool 
for obtaining Einstein equations but also for solving them. To this end we just note that the 
Bianchi identity takes in this formulation a very simple form: 

Da"^'^ a S^' = 0, (61) 

where \E'*-' is the matrix of self-dual components of the Weyl curvature tensor. In the case of 
the Schwarzschild solution the matrix \E'*-' is diagonal "^^^ = /9(r)diag(2, — 1, — 1), where /3(r) 



(57) 

The curvature = 
(58) 
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is a function of the radial coordinate. The Bianchi identity ( 16T1) then becomes a first order 
differential equation for /5, which immediately gives (3 ~ One then finds that one of the 

Einstein equations is an algebraic relation between g{r) and /3(r), so the problem of solving the 
system of Einstein equations reduces in this case to that of solving (I6T]) . This simple observation 
is relevant in situations other than Schwarzschild, and in some cases serves as an efficient tool 
for obtaining solutions. 

Finally, let us note that the described here method is closely related to that of Newman- 
Penrose [1], with the later being, loosely speaking, the two-component spinor version of the 
Plebanski one. In addition to working with self-dual quantities, and expressing all quantities 
in terms of spinors, the Newman- Penrose formalism introduces and works with a doubly-null 
tetrad, which is quite powerful, but at the same time makes the Lorentz-covariance of the theory 
not manifest. In contrast, Lorentz rotations of a tetrad in Plebanski formalism are described 
very simply as (complexified) SO (3) rotations. It is also considerably more easy to work with 
spacetime forms than with spinors. For these reasons Plebanski method may be an convenient 
middle ground between the relatively cumbersome tetrad method and a very powerful, but 
somewhat too relying on working in components Newman-Penrose formalism. 
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